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Abstract 



We present a determination of the gluon polarisation Ag/g'm the nucleon, based on the longitudinal 
double-spin asymmetry of DIS events with a pair of large transverse-momentum hadrons in the final 
state. The data were obtained by the COMPASS experiment at CERN using a 160 GeV/c polarised 
muon beam scattering off a polarised ^LiD target. The gluon polarisation is evaluated by a Neural 
Network approach for three intervals of the gluon momentum fraction Xg covering the range 0.04 < 
Xg < 0.27 . The values obtained at leading order in QCD do not show any significant dependence 
on Xg. Their average is Ag/g — 0.125 ±0.060 (stat.) ±0.063 (syst.) at Xg = 0.09 and a scale of 
ju2==3 (GeV/c)2. 
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1 Introduction 



The spin structure of the nucleon has been studied in polarised Deep Inelastic lepton-nucleon Scattering 
(DIS) for many years. The experimental observation by EMC that only a small fraction of the nucleon 
spin is carried by quark spins has strongly influenced more recent developments of spin physics. Several 
experiments were performed to confirm this result |2-8]. More measurements are in progress and/or in 
the data analysis phase: HERMES at DESY, STAR and PHENIX at RHIC, a number of experiments at 
JLAB, and COMR\SS at CERN. Several theoretical ideas were proposed |0] to explain this observation. 
In order to investigate the origin of the nucleon spin, it is essential to determine the spin fraction carried 
by gluons. Information about this quantity can be obtained indirectly from scaling violations in the 
structure function gi (see Refs. [s^ 10, JJ] and references therein) or from a direct measurement of the 
gluon polarisation (see Refs. 1 12 - lisll ). 



Leading order virtual photon absorption (LP) does not provide direct access to the gluon distribution 
since the virtual photon does not couple directly to the gluon. However, the observation of higher or- 
der processes opens a way to determine the gluon helicity distribution. Of particular interest is the 
Photon-Gluon Fusion (PGF) process shown together with leading-order photon absorption and QCD 
Compton scattering in Fig. [T] These processes are of first order in the strong coupling constant as, so 
their contributions to the DIS cross-section are comparable, but smaller than the virtual photon absorption 
contribution. 

The cleanest way to tag the PGF process is via open charm production, i.e. by selecting charmed mesons 
in the final state 1 18]. For this process the contribution from the leading order diagram is small because, 
in the COMPASS kinematic domain, the charm quark content in the nucleon is negligible. Due to the 
large mass of the charm quark, the contribution from fragmentation processes is also small. However, for 
the same reason, charm pair production in PGF is suppressed, so that the statistical precision on the gluon 
polarisation obtained in this way is limited. A way to overcome this limitation is to tag the PGF process 
leading to light quark pair production by detecting final state hadrons with large transverse momentum, 
Pt, with respect to the virtual photon direction. 

In the leading-order process, the hadron transverse momentum pj is due to the intrinsic transverse mo- 
mentum kr of quarks in the nucleon u% and to the fragmentation process, both resulting in small trans- 
verse momenta. A different situation occurs for QCDC and PGF processes, in which hadrons mainly 
acquire transverse momentum from the partons produced in the hard process. For this reason the re- 
quirement of observing two hadrons with large transverse momentum enhances the contribution of the 
PGF process in the selected sample 112011 . We present hereafter an analysis using this approach for the 
enhancement of PGF events in light quark production 12 iL 12211. 





a) b) c) 

Fig. 1: Feynman diagrams considered for y*N scattering: a) Leading order process (LP), b) gluon radia- 
tion (QCD Compton scattering), c) photon-gluon fusion (PGF). 
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3 DETERMINATION OFAG/G FROM MEASURED ASYMMETRIES 



2 Experimental set-up and data sample 

The experiment uses the naturally polarised muon beam at CERN. The experimental set-up consists 
of two major components: a polarised target and a magnetic spectrometer. A detailed description of 



the experiment can be found elsewhere [231] . A major upgrade of the COMPASS spectrometer was 
performed in 2005. For this analysis the most relevant improvement was a new target magnet which 
extended the angular acceptance. 

The average beam muon momentum is 160 GeV/c and the average beam polarisation is Ph = — 0.80± 
0.04 . The target consists of two cells in 2002-2004 and of three cells in 2006, located along the beam 
one after the other and filled with ^LiD. Lithium-6 can be regarded as a quasi-free deuteron and a helium- 
4 core. The average deuteron polarisation \Pf \ is about 0.5 and the average dilution factor of the target 
/ is 0.36. The latter is the ratio of the cross-section for all polarisable nucleons in the target material 
(deuterons) to that for all nucleons and includes radiative corrections. The relative uncertainties of |P,| 
and / are 5% and 2%, respectively. 

The data were collected during four years: 2002 to 2004 and in 2006. Selected events have an interaction 
vertex located in the target fiducial volume and contain both a beam muon and a scattered muon. The 
DIS region is selected by the requirement > I (GeV/c)^ and by a selection depending on the energy 
fraction y carried by the exchanged virtual photon, which leads to an invariant mass squared of the hadron 
system of > 5 (GeV/c)^. Events with y <0.l and with y > 0.9 are rejected because the former are 
more sensitive to experimental biases, while the latter are strongly affected by radiative effects. The 
above requirements define the inclusive sample. At least two additional charged hadrons associated 
with the vertex are required for the high-/??- sample. In the analysis the two hadrons with the highest 
Pt are selected and the following requirements are applied: > 0.7 GeV/c for the leading hadron, 
> 0.4 GeV/c for the sub-leading hadron, xp > for the Feynman variables of both hadrons and 
Zi +Z2 < 0.95, where zi 2 is the ratio of the hadron energy to the virtual photon energy. The cut on z 
removes events originating from exclusive processes. After all cuts, a sample of about 7.3 million events 
is used in the present analysis. 



3 Determination of Ag/g from measured asymmetries 

The longitudinal double-spin asymmetry for the production of two high-/7r hadrons in the DIS regime 
can be expressed as a function of the Bjorken scaling variable xbj. 

A^iixBj) = /?PGF«[g^ — (%) +RLFDA\''ixBj) +RQCDca2^''''A\''{xc) , (1) 

and all other variables are integrated over the experimental kinematic domain. The leading order (LO) 
inclusive asymmetry is given by the ratio of spin-dependent and spin-averaged quark distribution 
functions (PDFs), weighted by the squared quark electric charges; /?, is the fraction of process i and a'^^ 
the corresponding analysing power (i.e. the asymmetry of the partonic cross-section) {l^. The labels 
LP, QCDC and PGF refer to the processes presented in Fig. [T] The depolarisation factor D is the fraction 
of the muon beam polarisation transferred to the virtual photon and depends mainly on y. The variables 
Xbj, Xg and xc are the quark momentum fraction, the gluon momentum fraction in the PGF process and 
the quark momentum fraction in the QCDC process, respectively. Equation ([T|l is valid at LO in QCD 
assuming spin independent fragmentation. A possible spin dependence of fragmentation discussed in 
Ref. L25,1 can be neglected in the COMPASS kinematic region. 

The evaluation of Ag/g from the experimental asymmetry Af^ using Eq. ([T|) is possible only when the 
contributions from background processes (LP, QCDC) can be computed and subtracted. In this analysis, 
the fractions R, and the analysing powers a\^^ are extracted from Monte Carlo (MC). Therefore, the 
analysis requires a precise MC description of the data, so that /?/ and a^L calculated reliably. The 



5 



asymmetry can be evaluated from the spin-dependent and spin-averaged PDFs extracted from global 
fits or by using directly the measured inclusive lepton-nucleon asymmetry A^^ ■ In the present analysis 
we use the second option, which is less dependent on QCD analyses and related assumptions. As there 
are two unknowns in Eq. ([T]), A\'^{xBj) and A\'^{xc), the asymmetry A™^ has to be known for these two 
values of x and can be decomposed in a similar way as A^: 



(2) 



Combining Eqs. ^ and (O and neglecting small terms (note that the fractions /?pgf and /?qcdc are much 
smaller for the inclusive sample than for the high-;?/- sample), one obtains the following expression, 
which allows us to extract Ag/g: 



^PGFfl'LL^ — (%) 



+ 



R 



LP 



pincl 



incl.QCDC pincl 
«LL -^QCDC 



D 



pincl 
^LP 



pincl incl.PGpAg , 
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Here Eq. Q was used twice, once as given and once with the replacements Xg — )• x^, xc ^ and 

XBy Xc- 

Due to the fact that Ag/g is present in Eq. ([3]l at two different Xg values (denoted Xg and x'g), the extraction 
of Ag/g requires a new definition of the averaged Xg at which the result is obtained: 



Ai Xo — A2X' 

/ , ^ , where (4) 
M — M 

p incl p incl „QCDC 



1 _^PGFp incLPGFp "^PGF . j _ incl,PGF„ ^^PGF "LL 

^LP ^LP ^ 



Equation ^ relies on the assumption of a linear dependence of Ag/g upon Xg. The impact of the possible 
differences between Xg and x'g as well as between xc and x^ on the final Ag/g result is taken into account 
in the systematic uncertainty. 

The final relation between the gluon polarisation and A^l can be written as: 

Ag/g{x';) = "^"^Y^'^^ , with 



R 1 f a^^^^ ^incl.QCDC ^ \ 

a™"- = Al5;'(xB;)-^ +All (^c)^ ( -^j^Rqcdc - ^QCdc;^ j (5) 



incl,QCDC pincl „ QCDC 
^incl/^/ N ^LL ^QCDC «QCDC ^LL 



In the extraction of Ag/g we use a method similar to the one used in Ref. 1.261 . The target cells are 
labelled u, d for upstream and downstream. For 2006 the label u refers to the two outer cells and d to 
the central cell. The material in u and d cells is polarised in opposite directions. Spin orientations are 
reversed three times per day in 2002-2004 and once per day in 2006 by rotation of the target magnetic 
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4 MONTE CARLO OPTIMISATION AND NEURAL NETWORK TRAINING 



field by 180°. Data from before (u, d) and after such a rotation {u' , d') are combined in a so-called spin 
configuration, where nucleon spins in u and d' (d and u') have the same orientation. 

Data from different cells j = u,d,u' ,d' are combined so that beam flux, apparatus acceptance and spin- 
averaged cross-section cancel. The gluon polarisation is measured by solving the second order equation: 

PuPd' ^ (1 + (Ar)w + (A»)w Ag/g(x^-))(l + (A^?"-)>, + (A,,)>, Ag/g(x^-)) 

Pu'Pd (1 + (Ar>w + (A„')h. Ag/g(x-))(l + (A™-)„ + (Arf)„, Ag/g(x-)) ' ^ 

where pj is the sum of event weights w in sample j and {AJ"^)^ and {Aj)^ are weighted means of 
fPhPta^°" and fPi,Pt{X\ — A2), respectively. The weight w in the current analysis is defined as w = 
fPb{^\ — A2). In this way, Ag/g{x^g) is directly obtained, without going through the intermediate step of 
extracting the All(-^B7 ) asymmetry. 

In previous analyses of high-jo^ events iflil. \vh only mean values of Rj and a\j^/D for the three pro- 
cesses were used and the contribution of the leading process was suppressed by requiring the presence 
of two hadrons with high transverse momenta. Unfortunately, these requirements lead to a severe loss of 
statistics. In the present analysis, a Bayesian driven Neural Network (NN) approach for the extraction 
of Ag/g is used. It allows the use of loose pt cuts by dealing simultaneously with the three processes. 
The NN, trained on a MC sample, assigns to each event a probability to originate from one of these pro- 
cesses, which is then included in the weight w. Events more likely originating from processes other than 
PGF are kept with a small weight. For a given event, different NNs provide not only the probabilities 
to originate from a particular process but also the corresponding analysing powers and the momentum 
fractions xc and Xg. This approach makes optimal use of the data and avoids biases which may arise from 
correlations between analysing power and kinematic quantities used to evaluate the asymmetries. The 
statistical uncertainty of Ag/g is reduced by a factor of three comparing with the method used in [13]. 

4 Monte Carlo optimisation and Neural Network training 

In the present analysis the NN package from Ref. [27^ is used. Many results derived from a Neural 
Network approach strongly depend on the Monte Carlo sample on which the NN is trained. Thus, a good 
description of the experimental data by MC simulations is essential for the analysis. 

The LEPTO event generator 1I28I1 (version 6.5) is used to generate both an inclusive DIS sample and a 
sample which already contains at least two hi^-pT hadrons. The generated events were processed by 
the detector simulation program COMGEANT and reconstructed in the same way as real events by the 
reconstruction program CORAL. Finally, the same requirements are used in the analysis of real and MC 
events. 

Prior to the MC generator studies, an extensive effort was made to improve the detector simulation 
in COMGEANT. The MSTW08 PDF parametrisation ifjgll is used in the analysis as it gives reason- 
able agreement with measured in the COMPASS kinematic range [30J and is valid down to = 
1 (GeV/c)^. Also the function option from LEPTO is used, which improves data-to-MC agreement 
in the high-y region. Finally, a correction for radiative effects as described in Ref. [31] was introduced. 

The description of lepton variables was found to be satisfactory at this stage. For the hadron variables, 
the Parton Shower (PS) option in LEPTO had to be enabled to improve their description. However, this 
procedure introduces some inconsistency, since PS simulates higher order effects while the expression 
of Ag/g is derived at LO. The impact of this discordance will be taken into account in the evaluation of 
systematic uncertainties. In order to further improve the agreement with data for the hadron variables, 
some parameters describing the fragmentation process in LEPTO were tuned (high-;77 tuning in Table[T]). 
They correspond to the width of the gaussian pj distribution (PARI 21), the shape of the non-gaussian 
tail (PARI 23, PARI 24) and the symmetric Lund fragmentation function (PARI 41, PARI 42). 



7 



Table 1 : Default and tuned values of the LEPTO parameters describing the fragmentation process. 





PARJ21 


PARJ 23 


PARJ 24 


PARJ 41 


PARJ 42 


Default tuning 


0.36 


0.01 


2.0 


0.3 


0.58 


High- pt tuning 


0.34 


0.04 


2.8 


0.025 


0.075 




Cf {(GeV/cf) "bj y 

Fig. 2: Comparison between data (histogram) and MC simulations using high-/7r tuning (full squares) 
and default LEPTO tuning (open circles): distributions and Data/MC ratios for the lepton variables, Q^, 
XBj and y, normalised to the number of events. 

For the lepton variables the comparison of the high-p^ data sample to the MC sample is shown in Fig.|2] 
both for default LEPTO tuning and high-T^j- tuning. Figure [3] displays the corresponding comparison for 
the hadron variables (total and transverse momenta pi, pr^ of the leading and the sub-leading hadron 
P2, Pt2 and the hadron multiplicity. One observes that MC with high-/?7 tuning yields a satisfactory 
description of all distributions justifying its use to parametrise process fractions and analysing powers. 

Several NNs are used to parametrise all needed quantities. For a set of input parameters, the NN is trained 
to output the corresponding expectation value for a given quantity X. For the inclusive sample the input 
parameter space is spanned by xbj and Q^, while for the high-/77- sample the transverse and longitudinal 
momenta of the leading and sub-leading hadrons pxi , Pt2 > PLi > and pi^ are used in addition. 

An example of the quality of the NN parametrisation is given in Fig. |4l For the same MC sample it 
shows the probability for LP, QCDC and PGF events as a function of ^ p^ once as generated and once 
as obtained from the NN. The NN training was performed on a statistically independent MC sample. A 
good agreement is observed. While the LP probability reduces with increasing px (pr,, PT2 and Y^Pj), 
QCDC and PGF become the dominant contributions rising with similar strength. 



5 Systematic studies 

The main contribution to the systematic uncertainty comes from the dependence of the analysis on the 
MC. In total seven MC samples were prepared with different combinations of fragmentation parameters 
tuning (default LEPTO or high-pj), 'PS on' or 'PS off', different choices of the PDFs (MSTW08 or 
CTEQ5L is^] ) and Fl from LEPTO or from the R = Ol/Oj parametrisation of Ref. l33l] . In addition to 
what was already discussed, it is worth mentioning that for 'PS on' and 'PS off different so-called cut- 
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pJiGeWlc) 



2 4 6 8 10 12 14 16 1820 

Hadron multiplicity 



Fig. 3: Comparison between data (histogram) and MC simulations using high-T^j- tuning (full squares) 
and default LEPTO tuning (open circles): distributions and Data/MC ratios for the hadron variables, , 
PT2, Pi, P2 and the hadron multiplicity, normalised to the number of events. 
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Fig. 4: Values of /?lp, /?qcdc» ^pgf obtained from MC and from NN as functions of Y^pj (upper row), 
and their differences (bottom row). 



off schemes were used to prevent divergences in the cross-section calculations in LEPTO (see Ref. ll28ll ). 
These schemes and their parameters are quite important since their choice does not affect the data-to-MC 
comparison but changes the fraction of, e.g. PGF events. So, while keeping the default cut-off parameters 
proposed by the authors of Ref. tlstl . we tested various cut-off schemes. A small RMS value of 0.020 
was found for the Ag/g values obtained from these seven MC samples. However, it turned out that 
the asymmetry is very small, and so the above RMS may underestimate the systematic uncertainty 
related to MC. In order to avoid this, we consider in addition how the statistical uncertainty of Ag/g 
changes for various MC tunings. This leads to 5{Ag/g)Mc = 0.045 . 

The uncertainties of Ag/g due to the choice of the Af parametrisation and to the NN stability were found 
to be small, 5{Ag/g)^d = 0.015 and 5(Ag/g)NN = 0.010. The uncertainties of /, Pb, and P, have an 
even smaller impact on the final result: 5{Ag/g)f^Pi^ p^ = 0.004 . The HERMES results l34i] suggest that 
for heavier nuclei the dilution factor depends upon the transverse momentum of hadrons. Tests were 
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performed to check the dN jdpj^ dependence for the ^LiD target as compared to heUum, the medium in 
which the target material is immersed. No such dependence is observed. 

False asymmetries appear if the acceptance ratio of neighbouring target cells is different for the data taken 
before and after field reversal. They were searched for in a sample in which the event selection cuts were 
relaxed to to pj^-, > 0.35 GeV/c and > 0.7 (GeV/c)^. This leads to a large increase in statistics 
and allows for more precise studies of the spectrometer stability. No false asymmetries exceeding the 
statistical uncertainty were found. Taking the statistical uncertainty as limit for the false asymmetries 
one obtains 5(A^/g)faise = 0.019 . 

The two different values xq and xIq appearing in Eq. ([3]) were assumed to be equal. Two tests were done 
to check the systematic effect of this assumption. In the first one, a;^ = 1.6 •;cc was assumed, the value 
1.6 being an estimate taken from MC. In the second one, the NN parametrisation of xq was used with the 
previously obtained xq as input parameter instead of x^y. This leads to an uncertainty in Ag/g of 0.035 . 
Similar tests performed for x^ and changed x^^ by less than 0.01 . 

The expression used for the calculation of all assumes that the quarks are massless. This assumption is 
not valid for strange quarks. Tests were performed excluding kaons from the data sample, or making a 
parametrisation of the NN based on events with pions only. The final results are found to be stable within 
statistical fluctuations. 



Table 2: Summary of the contributions to the systematic uncertainty of Ag/^. 





Xg range 


[0.04,0.27] 


[0.04,0.12] 


[0.06,0.17] 


[0.11,0.27] 


MC simulation 


0.045 


0.077 


0.067 


0.129 


Inclusive asymmetry A\ 


0.015 


0.021 


0.014 


0.017 


NN parametrisation 


0.010 


0.010 


0.010 


0.010 




0.004 


0.007 


0.007 


0.010 


False asymmetries 


0.019 


0.023 


0.016 


0.012 


xq = x'(j in Eq. ^ 


0.035 


0.026 


0.039 


0.057 


Total systematic uncertainty 


0.063 


0.088 


0.081 


0.143 



The impact of resolved photon processes on the extracted value of Ag/g was studied using the RAPGAP 
generator [35']. It was found that events originating from resolved photons are expected to have very 
different kinematic distributions with respect to our standard high-/??- sample. It was checked whether 
adding an admixture of events originating from resolved photon processes would change the MC descrip- 
tion of the data. The results show that the contribution from resolved photons in our kinematic range is 
negligible. 

The contributions to the systematic uncertainty and their quadratic sum are presented in Table |2] They 
were also evaluated separately in the three x^ bins of Table [3] The total systematic uncertainty of the 
overall Ag/g result is obtained as 0.063, which is slightly larger than the statistical uncertainty. 



6 Results and conclusions 

The values of Ag/g provided by Eq. ^ were extracted for every spin configuration separateljQ in order 
to reduce systematic uncertainties. A correction for the probability of the deuteron to be in a D-wave 
state l,36il was applied. The mean values for each year of data taking are shown in Fig. |5] They are 



One configuration usually corresponds to 16h (2 days) of data taking in 2002-2004 (2006). 
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6 RES ULTS AND CONCL USIONS 



Table 3: Summary of the Ag/g results. 





Xg range 




[0.04,0.27] 


[0.04,0.12] 


[0.06,0.17] 


[0.11,0.27] 




0.09 


0.07 


0.10 


0.17 


Ag/g 


0.125 ±0.060 


0.147 ±0.091 


0.079 ±0.096 


0.185±0.165 



compatible within their statistical uncertainties and average to 



A^/g = 0.125 ±0.060 (stat.) ±0.063 (syst.) 



(7) 



at x^" = 0.09 and a scale of = 3 (GeV/c)^. 

The data cover the range 0.04 <Xg < 0.27 and were divided into three statistically independent subsam- 
ples in Xg as given by the NN. The con^elation between the generated Xg and the one obtained from the 
NN is about 62%. The results do not show any significant dependence of Ag/g on Xg (Table[3l). 

These results are compared with previous LO evaluations of Ag/g based on high-7?r hadron events in 
Fig. [6l The value taken from Ref. 1 17] is also derived from COMPASS data, however in the quasi-real 
photoproduction process instead of DIS. The hard scale and the range of gluon momentum are almost 
the same as in the present analysis and the two values of Ag/g are compatible within their statistical 
uncertainties. The Ag/g value obtained in the LO open-charm analysis [37] at a higher scale ^u^ = 
13 {GeV/cf is also shown. The SMC results from high-pr hadronpairs with > 1 (GeV/c)^ lEl] 



and the HERMES results from high-pr single hadrons using all II14I1 are compatible with the present 
results. 

The Ag/g curves shown in Fig. |6] are the results of global fits to spin asymmetries in inclusive and semi- 
inclusive DIS 1 3^ 3^. They were obtained at NLO in QCD and are thus not directly comparable with 
the LO result of the present analysis. It is however interesting to note that they all point to low values of 
A^/g for < 0.20. 




0.8 



2002 2003 2004 2006 



All 
data 



Fig. 5: Year by year Ag/g result and 
final average value. 



^ 0.6 
0.4 
0.2 


-0.2 
-0.4 
-0.6 
-0.8 



New COMPASS, high p^, Q^>1 (GeV/c)^, 2002-2006 
COMPASS, high p^, 0^<1 (GeV/c)^ 2002-2003 
COMPASS, open charm, 2002-2007 
SMC, high p^, 0^>1 (GeV/c)^ 
HERMES, high p , aii 




■ DSSV fit, n =3(GeV/c)'^ i 
LSS fit with AG>0, n^=2.5(GeV/c)^ 
LSS fit with AG changing sign, n^=2.5{GeV/c)^ 



10"' 



10-^ 



Fig. 6: Comparison of the final Ag/g with previous results 
(see text); the NLO curves are from Refs. |38, 39ll . 



A direct measurement of the gluon polarisation, extracted in the leading order approximation, was per- 
formed on all COMPASS data taken with a longitudinally polarised ^LiD target. The gluon polarisation 
Ag/g is extracted from a large sample of DIS events with > I (GeV/c)^ including a pair of high- 
Pt hadrons. A novel method using neural networks reduced the statistical uncertainty of the result and 
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allowed for the first time an evaluation of the gluon polarisation in three bins of the gluon momentum 
fraction Xg. 
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